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FOREWORD 
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of Dr. Robert Thomson, Head, Loads Control Section, Dynamics Branch. 

The technical monitor, Dr. Robert Hayduk, of this branch was the responsible 
government officer and assisted in some of the work. For the first two 
phases the co-principal investigators at VPI & SU were Dr. Robert Melosh 
and Dr. George Swift. Dr. Manohar Kamat, Mr. Ben Brenneman and the late 
Mr. Jon Dana provided technical support. For the third phase the co-principal 
investigators were Dr. Manohar Kamat and Dr. George Swift. Mr. Douglas Killian 
provided technical support on the second and third phases of the project. 

Dr. Manohar Kamat was the principal investigator for the final phases. 

Mr. Norman Knight, Jr., and Mr. Linh T. Duong provided technical support. 

The computer code CRASH previously developed by the Department of 
Transportation was the basis of the ACTION computer code which finally 
came to fruition only as a result of the efforts and dedication of several 
people, especially those of the graduate students. The author of this 
document is indebted to all of them for their contributions. The directions 
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and the continued encouragement of grant monitors Dr. Robert Thomson and 
Dr. Robert Hayduk of NASA Langley during the development of this program 
were as always invaluable. The hitherto taken-for-granted assistance 
of Ms. Barbara Durling of NASA Langley on several occasions during the 
development of the ACTION program certainly merits a great deal of ap- 
preciation. Finally, a word of gratitude is in order to Marlene Taylor, 
Fran Carter, B. J. Vickers and Jane Harrison for typing many versions 


of this document. 


ABSTRACT 


This document is a user's guide to the ACTION computer program. 

It defines the form and interpretation of input and output data. ACTION 
performs nonlinear transient response analysis of structures subjected 
to time varying loads, allowing for nonlinear, time independent material 
properties and large geometry changes. 
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SCOPE OF THE ACTION PROGRAM 


"ACTION" is an acronym for Analysis of Crash Transients in Inelastic 
and geometrically Nonlinear structures. The ACTION computer program uses 
the finite element method to determine the relationships between loads 
and deflections. 

Structures are modeled using a three-dimensional assemblage of 
line and sheet element types. Line elements include the rod, frame, and 
rigid types. Rod elements can only be loaded axially; frame elements can 
handle axial, bending, torsion and shear loads; rigid elements provide 
an infinitely stiff connection between two or more joints. The sheet 
element is a constant strain triangular membrane element capable of resisting 
both in-plane and out-of-plane loads through large changes in geometry. 

The program accepts both displacement and force boundary conditions. 
Concentrated forces may be applied at joints either as static loads or as 
time varying loads. Concentrated masses may also be located at joints. 
Displacement, velocity, acceleration and jerk (first derivative of 
acceleration) initial conditions are permitted. Any consistent system of 
units may be used for input data. The ACTION code includes a model of 
an impenetrable plane to simulate a flat, rigid ground plane, such as a 
concrete runway. Resistance to forward motion along this runway is provided 
by Coulomb friction. Any node of the model may become constrained to lie 
on the ground plane. Furthermore, any number of nodes may be attached to 
the ground plane at any one time. 
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ACTION predicts the response of a structure to static and/or dynamic 
loading conditions. The user selects an initial time step for transient 
problems, but the program determines subsequent time steps to satisfy error 
control criteria. The deformation model allows large changes to the 
initial geometry of the structure. The material model represents elastic- 
plastic response with strain hardening. It reflects cyclic loading including 
Bauschinger effects. The program also represents fracture of elements. 

Two algorithms are available for establishing equilibrium configurations: 
BFGS (Broden-Fletcher-Goldfarb-shanno) Variable Metric algorithm [1] and 
Powell's Conjugate Gradient algorithm [2], 

The BFGS algorithm permits static and dynamic analysis by minimizing the 
energy function. In the process it constructs data on the inverse hessian 
of the function and uses these data to accelerate the search process. Powell's 
conjugate gradient algorithm provides an alternate energy minimization- 
algorithm which, unlike BFGS, does not require storage for the upper triangular 
portion of the inverse hessian. Powell's conjugate method as implemented in 
ACTION has stricter convergence criteria than those of BFGS. Hence, difficult- 
ies in complete convergence may be experienced while using this algorithm. 

Output data consists of joint displacements, velocities, accelerations and 
jerks at the end of each time step. Stresses and strain rates are displayed 
at up to 16 locations throughout a cross-section. Elastic and dissipative 
energies are printed for each element. In addition, the total system energy, 
total elastic energy and total energy dissipation rate are printed. 

Data management features include capabilities to batch together static 
problems having identical geometry but different external loads, store data 
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on tape for subsequent restarts, restart from data stored on tape* Data 
may be stored on a restart tape by specifying a data dump for the end of a 
successful run. In addition, if calculation time approaches a user specified 
limit or if solution errors terminate a problem, data dumps are executed 
automatically. 

Figure 1 presents a general flow for the ACTION program. ACTION is 
mapped to handle 100 joints, 210 elements, five different types of material 
and 300 degrees of freedom. These limits may change, however, and the BLOCK 
DATA subroutine (not shown in Figure 1) should be consulted for a current 
specification of program limitations. In BLOCK DATA the variables NNODES, 
NELEM, NMATT and ND0F refer to the number of joints, elements, materials and 
the number of degrees of freedom for which the BFGS equation solver will be 
invoked is represented by NSOLVE in BLOCK DATA. 

This document cites the form and interpretation of input data required 
by ACTION and discusses the output data. The technical document [3] presents 
the mathematical basis of the code. 

This document is organized into three main sections. The next section 
describes the input data deck form and card entries. The second section 
illustrates the output and discusses its interpretation. The third section 
outlines a sample problem using some of the features of the code. 
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INPUT DATA 


The input data consists of a set of Data Blocks for each problem. 

The number of Data Blocks included in any particular problem is dependent 
upon the data required for that problem. The maximum number of Data 
Blocks that can be input per problem is ten. This section provides 
detailed instructions for preparing input. 

Input Data Deck 

Figure 2 shows how the Problem Data Decks form the Input Data Deck. 
Any number of problems can be batched together in one computer run. 

Two cards terminate the Input Data Deck. The next to the last card 
may be blank or contain any message the user desires. (It is convenient 
to use a message that indicates end of data for this computer run.) The 
last card is termed a Zero Card. This card may contain only zeros or 
blanks. (For ease of identification it is convenient to put a zero in 
all columns.) 

Problem Data Deck 

Figure 3 shows the format for each Problem Data Deck. The Problem 
Data Deck requires a Problem Title Card, three Control Parameter Cards, 
and at most ten Data Blocks. The Joint Coordinate and Material Data 
Blocks are required. The remaining Data Blocks are optional so long 
as they provide a sufficient definition of the problem. The Data Blocks 
included in the Problem Data Deck must appear in the order shown in 
Figure 3. 
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Figure 2 



Input Data Deck 
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GUESSED INITIAL DISPLACEMENTS I 


TERRAIN MODEL DATA BLOCK 


FORCING FUNCTION DATA BLXK 


OPTIONAL 

DATA 

BLOCKS* 


FUNCTION TABLE DATA BLOCK 


INITIAL CONDITION DATA BLOCK 


MASS DATA BLXK 


MEMBRANE ELEMENT DATA BLXK 


LINE ELEMENT DATA BLXK 


MATERIAL DATA BLXK 


JOINT CXRDINATE DATA BLXK 


CONTROL PARAMETER CARDS (1) 


PROBLEM TITLE CARD 


DATA BLXKS SELECTED 
MUST BE SUFFICIENT TO 
DEFINE A PROBLEM 


REQUIRED 

INPUT 

DATA 


Figure 3. Problem Data Deck 




Figure 4. Data Block Organization 
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1 . 


Problem Title Card : One card for each problem. The card contains 


any alphameric message desired by the user. The format is given 
in Table 1. 

2. Control Parameter Cards : Three cards for each problem. The first 

card contains information used in selecting program options. The 
second card contains information used to control the step-by-step 
numerical integration in time. The third card contains information 
about the 3 and y parameters of the Newmark-Beta method as well as solu 
tion and function accuracy limits. The formats are detailed in Table 2 

3. Data Blocks : Each Data Block contains information for describing 

a particular feature of the problem at hand. Figure 4 indicates the 
anatomy of each Data Block. Each Data Block begins with a Name 
Card and ends with a Zero Card. The first four columns of a Name 
Card must duplicate exactly the first four letters of the Data 
Block name as it appears in Figure 3. The remaining columns may 
contain any message desired by the user. The Zero Card may contain 
only zeros or blanks. (For ease of identification it is convenient 
to put a zero in all columns.) Data Cards peculiar to each Data 
Block are placed between the Name Card and Zero Card. If all values 
for an entire Data Block are zero, the entire Data Block including 
the Name Card and the Zero Card can be omitted from the Problem Data 
Deck. 

Tables 3 through 12 define and interpret the input data required 
for the various Data Blocks. These tables have been made self-contained 
through the use of footnotes and figures. In these tables !f b M is used 
to indicate a blank (unpunched) card field. Card formats are listed 
using FORTRAN notation. 
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Table 1 


PROBLEM TITLE CARD 


Card 

Cols 

Format 

Item 

1 

i 


" b " 


2-80 

A 

Problem Title 1 


: An alphameric problem title prescribed by the user. This title is used 
as the heading of the problem printout. 
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Table 2 


CONTROL PARAMETER CARDS 


Card 

Cols 

Format 

Item 

i 

1-5 

15 

IRUN, source of input data 
1, read input from cards 
-2, pertains to tape 10 1 


6-10 

15 

I INK, print level for minimization progress 

0 , minimum 

1, standard 

2, detailed 


11-15 

15 

ISINK, stress print flag 

0, no stresses 

1, stresses 

2, stresses, including detailed printout 
for membrane elements 


16-20 

15 

ITAPE, receptacle for storing restart data 
0, no restart data need be saved 
-1, pertains to tape 10 2 


21-25 


Not used 


26-30 

15 

ILIN, analysis sophistication 

0, infinitesimal deformations (linearity) 

1, infinitesimal joint motions and finite 
element distortions (small strains) 

2, finite joint motions and element 
distortions (small strains) 

3, finite joint motions and element 
distortions (moderately large strains) 3 


31-35 

15 

IC0NT, batching of static problems with 
identical geometry but different external 
loads 4 

0, batching feature not desired 

1, batching desired 

- 

36-40 

15 

MLIN, maximum number of linear minimizations 
(Defaults to three times the number of degrees 
of freedom if entered as zero) 

- 

41-45 


Not used 
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Table 2 (Continued) 


CONTROL PARAMETER CARDS 


Card 

Cols 

Format 

Item 


46-50 

15 

I ORDER (Defaults to 2) 

1, first order method (to be used when 
IS0LVE = 5) 

2, second order method (to be used when 
I SOLVE = 7) 


51-55 

15 

IS0LVE, equation solver option (Defaults to 7) 
7, BFGS method used unless dynamic memory 
allocation overflows available storage - 
then Powell’s method automatically 

invoked ^ 

5, Powell’s conjugate gradient method 
used entirely 


56-60 


Not used 


61-65 

15 

IPRT, Interval of time steps printouts 
(Defaults to 1 - every time step is printed. 
With IPRT-99999 only the results for the 
first and the last time are printed.) 


66-75 


Not used 


76-80 

15 

ISPEC 

0, No non-zero prescribed displ. 

1, Some non-zero prescribed displ. 


^Restart data was written onto tape 10 during a previous run and is to be used 
as input for the current run. For problems restarted from tape, the Problem 
Data Deck consists of only the Problem Title Card and the three Control 
Parameter Cards. 

2 In order to store data on tape 10, a scratch work space, tape 12, is 
required. Tape 12 is also required as a scratch work space to handle 
increased dimension of the inverse hessian approximation when frame elements 
yield. This is due to an additional axial degree of freedom used for the 
definition of the deformation characteristics of a frame element that has 
yielded. 

^This option should be invoked only rarely at the discretion of the user, 
since it does not account for large strains in a consistent manner. The 
only thing that is achieved by invoking this option is the use of nonlinear 
strain displacement relations in the corotational co-ordinate system. No 
distinction is made between pseudo and Cauchy stresses nor are changes in 
density accounted for [3]. 
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Table 2 (Continued) 


CONTROL PARAMETER CARDS 


Card 

Cols 

Format 

Item 

2 

1-10 

E10.0 

TBEGIN, time to start integration 


11-20 

E10.0 

HALT, time to stop integration 


21-30 

E10.0 

TSTEP, starting time step size 


31-40 

110 

NTIMES, maximum number of time steps or 
load steps to take 5 


41-50 

E10.0 

BTERR, biggest tolerable error (a larger 
error causes a decrease in time step) 5 


51-60 

E10.0 

STERR, smallest tolerable error (a smaller 
error causes an increase in time step) 5 


61-70 

E10.0 

1 

TMIN, minimum time step permitted 


71-80 

E10.0 

ELAPS, maximum running time per job as 
measured by CPU clock, minutes (Defaults 
to 100 CPU minutes) 


^A "batched" problem is entered as a "pseudo-transient" problem with zero 
displacement initial conditions. The only distinction in execution is that 
displacements are set equal to zero after each load step. Loading conditions 
are obtained from Function Tables. See footnote 5 for an explanation of 
"pseudo-transient" problems. 

5 A transient problem for which no dynamic response is desired is called a 
"pseudo-transient" problem. This may be used for the "static analysis" 
of a structure subjected to incremental changes in load. Loads are input 
using Function Tables, but force initial conditions must duplicate forces 
from Function Tables for the TBEGIN time. Pseudo-transient problems 
require that: 

a) NTIMES be equal to the number of loads in addition to the initial 
load for which a static analysis is required. 

b) No masses be input. 

c) Function Table time steps must be constant and equal to TSTEP. 

d) BTERR must be equal to zero to make integration time steps 
consistent with Function Table time steps. 

^The selection of suitable error bounds requires judgment on the part of 
the user. Too large a time step results in truncation error (displacement 
response is assumed to be a third order function of time) ; too small a 
time step reduces the accuracy of the computer representation of the 
difference in displacements over the length of the time step. 
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Table 2 (Continued) 


CONTROL PARAMETER CARDS 


Card 

Cols 

Format 

Item 

3 

i-io i 

E10.0 

BETA, Newmark-Beta Integration Parameter 
(Defaults to 0.276 [4]) 


11-20 

E10.0 

GAMMA, Newmark-Beta Integration Parameter 
(Defaults to 0.55 [4]) 


21-30 

E10.0 

DISACC, Desired Solution Accuracy 
(Defaults to l.E-08) 


31-40 

E10.0 

FUNACC, Estimated Function Accuracy 
(Defaults to l.E-13) 


41-50 

E10.0 

TSTEP1, Restart with a step size which 
is required to be different from the 
automatic step size determined by the 
analysis 7 


7 This feature can be 
analysis at the end 


used only when restarting from a successfully completed 
of a time step. 
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Table 3 


JOINT COORDINATE DATA BLOCK 1 
(required) 


Card 

Cols 

Format 

Item 

■ 

i 

1-80 

A 

"JOIN ..." 


1-5 

15 

Joint number 


6-10 


Not read 


11-40 

3E10. 

The X, Y, and Z coordinates of the joint 


41-44 


Not read 

u 

a 

45 

11 

X-displacement 1 

•H 



fl 

O 



1 

i“> 

46 

11 

Y-displacement 1 Restraint codes 

Xl 

CJ 



S 0 = unrestrained 

d 

<u 

47 

11 

Z-displacement I 1 = restrained 

u 



\ 

o 

48 

11 

X-rotation / 

'd 



| 

u 

aj 

49 

11 

Y-rotation | 

o 




GJ 

d 

50 

11 

Z-rotation ■ 

o 

51-80 


Not read 

Last 

1-80 


"00000 ..." or "bbbbb ..." 


lr The joint coordinate data block consists of a data block name card, one 
card for each joint and a zero card to indicate the end of the data block. 
The joints must be numbered sequentially beginning at one; no number may 
be skipped or omitted; the total number of joints equals the largest 
joint number. Joint data cards must be in sequential order. 
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Table 4 


MATERIAL DATA BLOCK 1 


Card 

Cols 

Format 

Item 

i 

1-80 

A 

"MATE ..." 

1 

1-5 

1.5 

Material number 


6-10 


Not read 


11-80 

A 

i 

Material name (User f s Option) 

2 

1-80 

8E10.0 

a 1 , e , a , e„, a„, e_, a,, e, / stress-strain 




V curve 


1-80 

8E10.0 

a^, e^, o^, £y , Og, £g l coordinates 2 

4 

1-20 

2E10.0 

V V 


21-30 

E10.0 

Elastic shear modulus 


31-40 


Not read 


41-50 

E10.0 

Mass Density^ 


51-80 


Not read 

Last 

1-80 


"00000 ..." or "bbbbb ..." 


1 The materials data block consists of a data block name card, then four 
cards for each type of material to be specified, and a zero card to 
indicate the end of the data block. 

2 For each material, the uniaxial stress-strain curve is modeled by eight 
straight-line segments, as shown in Figure 5a. Four segments must be used 
for the tensile side and four for the compressive side; fewer than eight 
segments are not allowed . 

For eg ± e <_ the material is treated as elastic. Yielding occurs 
at £2 tension and eg in compression. The elastic modulus in tension may 
not be different from the elastic modulus in compression. (The elastic 
moduli are the slopes of the elastic segments of the stress-strain curve.) 
The stress-strain curve must be single-valued in strain. This requires 
that the strain continuously increase with distance from left to right along 
the curve. However, it is permissible to specify an elastic perfectly 
plastic material. The specification of a rigid work-hardening material is 
not possible since the slope of the elastic segments of the curve cannot 
exceed the largest number than can be accommodated by the computer in 
question. The slope of the first plastic segment of the curve in tension 
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and compression must be less than the corresponding elastic tension or 
compression slope. Otherwise, there is no restriction on the slopes of 
the stress-strain curve, except the stress must be uniquely defined 
for any given strain. 

For unloading from the plastic segments of the stress-strain curve, the 
material response is elastic as described by the elastic segments of the 
curve. This is shown by the dotted lines. Note that the segments describing 
initial elastic response are translated to the plastic unloading point and 
used to describe the elastic unloading. Thus, strain hardening and 
Bauschinger effects can be represented. 

It is noted that the point (o]_, e^) of the material stress-strain curve 
may be displaced from the origin to simulate special conditions of initial 
stress or strain that may not be simulated by the initial conditions 'data. 

For inadmissible, indeterminate stress-strain conditions, see Figures 5b and 
5c. In these figures points on certain portions of the curves are not uniquely 
defined. 

^Nonzero density provides an automatic calculation of the lumped masses at 
nodes. If the user chooses to specify his own lumped mass distribution, 
the mass density should be set to zero. Lumped masses in addition to 
those automatically calculated by the program as a result of the nonzero 
density specification, may also be specified by using the MASS DATA BLOCK 
feature to be described later. 


a 



Figure 5a. Stress-strain Curve Representation 
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STRAIN 


Figure 5b. Indeterminate Stress-Strain Conditions 



STRAIN 


Figure 5c. Indeterminate Stress-Strain Condition 
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For Each 


Table 5 


LINE ELEMENT DATA BLOCK 1 
(may be omitted) 


Card Cols Format 


Item 


1 

1-80 

A 

For Each 

9 

Member 
oj ro h* 



Last 

1-80 



"LINE . . . 


Line Element Parameter Card 3 
Cross Section Data 11 
Cross Section Data 1 * 


"00000 . . . " or "bbbbb . . . 


^he line element data block consists of a data block name card, three data 
cards for each member, and a zero card to indicate the end of the data block. 

^Element data can be put in any order, but cards for each element must be in 
order. 

3 Refer to Table 5.1. 

4 The information on the two cross section data cards depends on the type of 
member . 

Rod Member - Refer to Table 5.2. 

Frame Member - General Cross Section - Refer to Table 5.3. 

Frame Member - IE Section - Refer to Table 5.4. 

Frame Member - Box Section - Refer to Table 5.5. 

Frame Member - Circular Tube Section - Refer to Table 5.6. 

Frame Member - Elliptical Section - Refer to Table 5.7. 

Frame Member - Solid Rectangular Cross Section - Refer to Table 5.8. 

Rigid Member - Refer to Table 5.9. 
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Table 5.1 


Cols 

1-5 

6-10 

11-15 

16-20 


21-25 


26-30 

31-35 

36-40 

41-45 


LINE ELEMENT PARAMETER CARD 


Format 


Item 


15 

15 

15 

A4,1X 


15 


15 


15 


15 

15 


Member number . 

The joint number for one end of the member, henceforth 
the P-end , or primary node for a STIF member. 

The joint number for the other end of the member, 
henceforth the Q-end (or total number of secondary 
nodes for a STIF member) . 


Type of member. 
Alphameric Code ' 


Description 


ROD 

Rod Member 

GLOB 

Frame Member of General Cross Section 

IE 

Frame Member, IE Section 3 

BOX 

Frame Member, Box Section 4 ' 

TUBE 

Frame Member, Circular Tube Section 5 

ELIP 

Frame Member, Elliptical Section 6 

SORE 

Frame Member, Solid Rectangular 


Section 7 

STIF 

Rigid Member 8 


Member status for rod and frame members. 

1 - Unimodular materials in elastic range 

2 - Initially unimodular and elastic, but may change 

to plastic 

3 - Requiring numerical integration over the cross 

section 

If the element is assumed to lie along a plane of 
symmetry the member status is specified with a negative 
sign, i.e., -1, -2 or -3 as required. 

Material number. 

Not required for rigid members. 

Code parameter for rod and frame members. 

Negative, Compression only (gap or snubber) 

0, Normal Member 

Positive, Tension only (guy wire) 9 
Not used 

Shear deflection for frame members 10 

0, Shear deflection of frame member neglected 

1, Shear deflection of frame member permitted 
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Table 5.1 (Continued) 


LINE ELEMENT PARAMETER CARD 11 


Cols 

Format 

Item 

46-50 


Not read 

51-60 

E10.0 

Angle <j) (in degrees) which describes the roll 

X 

* 


orientation of frame members. See Figure 6. 

61-80 


Not read 


1 Refer to Figure 6. 

Alphameric codes can be punched anywhere in the field but embedded blanks 
are not allowed. 

3 Refer to Figure 7. 

4 Ref er to Figure 8. 

5 Refer to Figure 9. 

6 Refer to Figure 10. 

7 Refer to Figure 11. 

8 A rigid member is an infinitely stiff member. It can be used to represent 
rigid links in the model. 

9 Tension only is not presently an option, but it can be obtained in rod 
members by using a zero compressive load. 

^ Shear deflection for nonlinear members should be requested only when 
needed, since shear calculations may increase computer run time 
significantly. 

1 ! Any data not required of a particular member type should be entered as 
zero or left blank. 
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Table 5.2 


ROD MEMBER 


Card 

Cols 

Format 

Item 

2 

1-10 


Not read, it is convenient to place the 
member type in Cols 1-10. 


11-20 

E10.0 

Cross section area. 


21-30 

E10.0 

Cross section moment of inertia for 
buckling (usually the minimum transverse 
inertia) . 


31-80 

5E10.0 

Not used. 

3 

1-10 


Not read, it is convenient to place the 
member type in Cols 1-10. 


11-80 

7E10.0 

Not used, but this card must be present. 
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Table 5.3 


FRAME MEMBER OF GENERAL CROSS SECTION 1 


Card 

Cols 

Format 

Item 

2 

1-10 


Not read, it is convenient to place the member 




type in Cols 1-10. 


11-80 

7E10.0 

Area, I B > I BC , TJ, Sg, S ( ,. 1 

3 

1-10 


Not read, it is convenient to place the member 




type in Cols 1-10. 


11-30 

2E10.0 

v V 


31-80 

5E10.0 

Not used. 


1 Cross section parameters are described with respect to axes B and C which 
are situated in the plane of the cross section with their origin at the 
centroid. These axes need not be principal axes of the cross section. 

AREA = Cross section area 

I = Area moment of inertia about B axis 
B 

I = Area moment of inertia about C axis 
L/ 

I = Area cross product of inertia 
BC 

TJ = Cross section geometrical torsion constant 

S = Distance to shear center along B axis 
B 

S p = Distance to shear center along C axis 
= Shear deflection constant, B direction 
K = Shear deflection constant, C direction 
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Table 5.4 


FRAME MEMBER-IE SECTION 1 


Card 

Cols 

Format 

Item 

2 

1-10 

11-80 

7E10.0 

Not read, it is convenient to place the member 
type in Cols 1-10. 

d-j^ , t.^, d^ 9 ^ 3 * ^4 

3 

1-10 

11-80 

7E10.0 

Not read, it is convenient to place the member 
type in Cols 1-10. 

V d 5’ t 5’ d 6 ’ t 6* d 7 ' t 7 2 

2 Ref er 

to Figure 

7. 



2 For the IE-Section any or all of the five flange pieces may be omitted. 
The two web pieces — t^, and d^ s t^_ — must be non-zero. 
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c, centroid 

s, shear center 

Equally spaced stress points 
in each flange or' web piece. 


STRESS POINT 



Table 5.5 


FRAME MEMBER-BOX SECTION 1 


Card 

Cols 

Format 

Item 

2 

1-10 


Not read, it is convenient to place the member 
type in Cols 1-10. 


11-50 

4E10.0 

d l» t v d 2’ t Z 


51-80 

3E10.0 

Not used. 

3 

1-10 


Not read, it is convenient to place the member 
type in Cols 1-10. 


11-80 

7E10.0 

Not used, but this card must be present. 


1 Refer to Figure 8. 
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Table 5.6 


FRAME MEMBER-CIRCULAR TUBE SECTION 1 
Card Cols Format Item 


2 

1-10 


Not read, it is convenient to place the member 




type in Cols 1-10. 


11-30 

2E10.0 

d, t 


31-80 

5E10.0 

Not used. 

3 

1-10 


Not read, it is convenient to place the member 




type in Cols 1-10. 


11-80 

7E10.0 

Not used, but this card must be present. 


1 Refer to Figure 9. 
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Table 5.7 


FRAME MEMBER-ELLIPTICAL SECTION 1 


Card 

Cols 

Format 

Item 

2 

1-10 


Not read, it is convenient to place the member 
type in Cols 1-10. 


11-50 

3E10.0 

d i* d 2 


51-80 

4E10.0 

Not used. 

3 

1-10 


Not read, it is convenient to place the member 
type in Cols. 1-10. 


11-80 

7E10.0 

Not used, but this card must be present. 


1 Refer to Figure 10. 

2 If di and d .2 are very nearly equal the section should be approximated by a 
tube of average diameter. 
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Table 5.8 


FRAME MEMBER-SOLID RECTANGULAR SECTION 1 


Card 

Cols 

Format 

Item 

2 

1-10 

11-50 

51-80 

2E10.0 

4E10.0 

Not read, it is convenient to place the member 
type in Cols 1-10. 

d i» fc i 
Not used. 

3 

1-10 


Not read, it is convenient to place the member 
type in Cols. 1-10. 


11-80 

7E10.0 

Not used, but this card must be present. 


1 Refer to Figure 11. 
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Figure 11. Geometry of Solid Rectangular Cross Section 
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Table 5.9 


RIGID MEMBER 


Card 

Cols 

Format 

Item 

2 

1-10 

11-80 

7E10.0 

Not read, it is convenient to place the member 
type in Cols 1-10. 

Secondary node numbers 1 

3 

1-10 


Not read, it is convenient to place the member 
type in Cols 1-10. 


11-80 

7E10.0 

Not used, but this card must be present. 


^The only data required for the rigid members are the numbers of the nodes 
connected together by the rigid member. Card 1 contains the number of 
the primary node as indicated previously and card 2 contains the numbers 
of the secondary nodes. The primary node is the only one to which degrees 
of freedom for rigid body motion are assigned. The choice as to which of 
the nodes is the primary one is arbitrary. 
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Table 6 


MEMBRANE ELEMENT DATA BLOCK 1 


Card 

Cols 

Format 

Item 

i 

1-80 

A 

"MEMB ..." 


1-5 

15 

Member number 


6-10 

15 

P -joint number | 


11-15 

15 

Q-joint number > P, Q, R specified in a 




V counterclockwise manner 2 


16-20 

15 

R-joint number I 

u 



' 

o 




X 

B 

21-25 

A4,1X 

Alphameric member code - MBRN 

CD 




B 




X 

26-30 

15 

Member status 

o 

CT3 



1, Elastic only 

CD 



2, Initially elastic, but may go inealstic 

U 

o 



3, Inelastic 

4-1 




X 

u 

31-35 

15 

Material number 

d 




o 




0) 

36-40 

15 

Not read 

G 




o 

41-50 

F10.0 

Thickness of the membrane element 


51-80 


Not read 

Last 

1-80 


"00000 ..." or "bbbbb ..." 


1 The membrane element data block consists of a data block name card, one 
card for each member and a zero card to indicate the end of the data 
block. 

9 

The three stress and strain components are calculated at the centroid of 
the element with reference to the local element coordinate system described 
in Ref. [3] . 
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Table 7 


Card 


Cols 


MASS DATA BLOCK 1 
(may be omitted) 

Format 


Item 


1-80 

A 

1-5 

15 

6-10 


11-70 

6E10.0 

71-80 


1-80 



"MASS 


CO 

CO 

cd 

S 

A 

u 

■H 

£ 

U 

a 

•H 

O 


-C 

a 

td 

<u 

o 

4-1 

"d 

u 

cd 

o 

a) 

a 

o 


Joint number 
Not read 

M^, Mass parameter for motion in the X-direction 

My* Mass parameter for motion in the Y-direction 

M^, Mass parameter for motion in the Z-direction 

1^, Mass moment of inertia for rotation about 

X-axis 

I , Mass moment of inertia for rotation about 
Y-axis 

1^, Mass moment of inertia for rotation about 
Z-axis 


Not read 


Last 


"00000 . 


. " or "bbbbb . 


*The mass data block consists of a data block name card, one card for each 
joint with non-zero lumped masses or inertias, and a zero card to indicate 
the end of the data block. The program does not treat products of inertia. 


2 If M x is negative, M^, My and are set equal to | M^. | . In this case My 
and are ignored, regardless of their values. 
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Table 8 


INITIAL CONDITION DATA BLOCK 1 
(may be omitted) 


Card 

Cols 

Format 

Item 

i 

1-80 

A 

"INIT ..." 


1-5 

15 

Joint number 

CO 

6-10 

15 

Type of initial condition 2 

u 

a 



1. Displacement 

•H 

o 



2. Velocity 




3. Acceleration 

u 

o 



4. "Jerk" 

M-l 



5. Force 





o 




u 

•H 

11-40 

3E10. 

Translation component, X, Y, Z-direction 

3 

cr 




a 

u 

41-70 

3E10. 

Rotation component, X, Y, Z-axes 

CO 

< 

71-80 


Not read 

Last 

1-80 


"00000 ..." or "bbbbb ..." 


J The initial condition data block consists of a data block name card, one 
card for each joint with non-zero initial conditions, and a zero card to 
indicate the end of the data block. 


2 Multiple initial conditions may be treated by using multiple cards with the 
same joint number. The program will superimpose all conditions given. 
However, care must be exercised to avoid inconsistent conditions. That is, 
an initial force implies an initial displacement and, if the force is 
suddenly removed, an initial acceleration. Thus, if an initial force is 
specified, the program calculates the associated displacements and 
accelerations. Normally other displacement or acceleration values should 
not be specified by the user. If initial displacements are non-zero, the 
program calculates the associated initial stresses. 

The initial forces introduced here may be regarded as the forces 
required to produce initial displacements. The option of specifying the 
forces, rather than the displacements, relieves the user of the task of 
calculating these displacements. These are static forces applied prior 
to the start of the integration to produce static, initial displacements. 
When the integration starts, these forces are removed. If the problem 
requires that the force be maintained after the start of integration, 
the forcing function feature should be used to maintain the desired force 
after time integration begins. 

The numerical integration procedure is such that initial accelerations 
and "jerks" (first derivative of acceleration) should be specified if 
applicable. This is contrary to the familiar analytical integration of 
equations of motion where only initial displacements and velocities are 
required. 


38 








Table 9 


FUNCTION TABLE DATA BLOCK 1 
(may be omitted) 


Card Cols Format Item 


1 

1-80 

A 

"FUNC ..." 

1 

1-5 

15 

Function number 


5-10 

15 

Number of time points 


11-80 

A 

Function name (user’s option) 


1-10 


Not read 2 

QJ 

U 

11-70 

71-80 

6E10.0 

tf> F-^ 5 , • • • 

Not read 

2 

u 4 

<D 

U 

1-10 


Not read 2 

t n 
< 

11-70 

6E10.0 

• 

* 

rt 

2 

2 


71-80 


Not read 

Last 

1-80 


"00000 ..." or "bbbbb ..." 


*The function table data block consists of a data block name card, data 
cards for each function to be defined, and a zero card to indicate the 
end of the data block. 

The function tables are used to define reference time-functions for 
describing time-varying joint forcing functions. 

Each time function is defined by 
straight-line segments connecting pre- 
scribed time points, as sketched to the 
right. There is no explicit limit on 
the number of time points. There must 
be at least one time point. 

There is an implicit limit on the 
number of time points. The function 
tables data is stored dynamically in a 
data array along with other program 
data. When the array size is exceeded 
an error message is printed. 
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The function must define a continuous increase in time; "backward" slopes 
are not permitted. Vertical slopes are acceptable. For vertical slopes 
the average value of the two end points is used. 

The functions must be numbered sequentially beginning at one end; no 
number may be skipped or omitted; the total number of functions equals 
the largest function number. The function data cards must be ordered 
sequentially by function number. 

2 It is convenient to place the function number in Cols 1-5. 
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Table 10 


FORCING FUNCTION DATA BLOCK 1 
(may be omitted) 


Card Cols Format Item 


1 

1-80 

A 

"FORC ..." 


1-5 

15 

Joint number 


6-10 

15 

Degree-of-freedom parameter 



i 

1, Force, x-direction 

2, Force, Y-direction 

3, Force, Z-direction 




4, Moment, x-direction 

o 

*H 



5, Moment, Y-direction 

U 

o 

G 



6, Moment, z-direction 

3 

4-1 

11-15 

15 

Reference function parameter 

toO 

a 



Positive, Function table number 

*H 

a 



-1, Constant function of magnitude A3 

u 

o 



-2, Sinusoidal function 

4-1 



-3, Cosinusoidal function 

rC 

a 




ej 

Q> 

16-20 


Not read 

u 

o 

21-60 

4E10.0 

Al, A2 , A3, A4 




Al, Scaling parameter (multiplication) 

u 

cd 



A2, Time shift parameter 

a 



A3, Magnitude shift parameter (addition) 

<u 

a 



A4 , Time period for sinusoidal or 

o 



cosinusoidal functions, not used 
otherwise. 


61-80 


Not read 

Last 

1-80 


"00000 ..." or "bbbbb ..." 


*The forcing function data block consists of a data block name card, one 
card for each forcing function to be prescribed, and a zero card to 
indicate the end of the data block. This data block defines time-varying 
forcing functions applied at the joints. The time variation is described 
by a reference function. The reference function may be obtained from the 
function tables, may be constant, or may be sinusoidal or cosinusoidal. 
The reference function may be scaled by a multiplication factor (Al) , 
shifted in time by a time factor (A2) or shifted in magnitude by an 
addition factor (A3). Refer to the Fig. 12 on the next page. Multiple 
forcing functions may be superimposed at a single degree-of-f reedom. 
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Table 11 


TERRAIN MODEL DATA BLOCK 1 
(may be omitted) 


Card Cols Format Item 


1 

1-80 

A 

"TERR ..." 

2 

1-10 

E10.0 

M Y" coordinate of ground plane 


11-20 

E10.0 

Coefficient of friction of ground plane 


21-80 


Not read 

As re- 
quired 

1-5 

15 

Number of the joint initially constrained 
to lie on the ground plane 




Not read 

Last 

1-80 


"00000 ..." or "bbbb . . . " 


1 The terrain model data block consists of a data block name card, a card 
describing the ground plane, a card for each joint initially on the ground 
plane, and a zero card. 

This data block is used to establish a plane which the structure is 
not allowed to penetrate. The structure may slide along the plane. The 
plane is assumed to be the X-Z plane of the global reference system 
located by the Y coordinate as input. 
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Table 12 


GUESSED INITIAL DISPLACEMENTS DATA BLOCK 1 
(may be omitted) 


Card Cols Format Item 


1 

1-80 

A 

"GUES ..." 

As 

required 2 

1-80 

4E20.0 

Initial guess for displacement vector, 
be read in order of degrees of freedom 
must include all degrees of freedom- 3 

Must 

and 

Last 

1-80 


"00000 ..." or "bbbbb ..." 


*The guessed initial displacements data block consists of a data block name 
card, a set of data cards specifying an initial guess for displacements 
and a zero card to indicate the end of the data block. This data block 
is used to provide the minimization process with an estimate of the 
displacement solution vector. Minimization begins from this initial guess. 

2 The number of data cards required is equal to the number of degrees of 
freedom divided by four. 

3 This feature can be used only if initial conditions are specified, that 
is if an Initial Condition Data block is present in the input. 
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OUTPUT DATA 


Output data includes both printed data and data stored on 
peripheral storage units. Printed output includes an interpretation 
of problem input, a record of the time response of the structure, and 
error messages. Data stored on peripheral storage units may be used to 
restart the analysis or plot the structural response. This section 
presents a description and interpretation of the output data. 

Printed Output 

Printed output consists of a definition of the problem as obtained 
from input data, a description of equilibrium positions for each time 
or load step, including a record of minimization progress, and diagnos- 
tic messages printed when anomalies arise during execution. 

The amount of printed output is controlled by the IINK and ISINK 
parameters on the first Control Parameter Card in Table 2. IINK selects 
the minimization print level and ISINK controls stress printout. Some 
printout, however, cannot be controlled by the user. Joint displace- 
ments, velocities, accelerations and jerks, energy summaries and error 
measures are always printed at the end of first and last time or load 
step and at specified intervals. 

The following pages present figures illustrating typical printed 
output. These figures do not represent a sample problem, but include 
output from several problems in order to illustrate as many features of 
the output as possible. The computer page representations in these 
figures are reduced reproductions of the actual computer printout. 

There are a maximum of 120 characters per line of printout. Typewritten 
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notes do not appear on the computer printout but have been added to 
provide a self-contained interpretation of output data. 

Actual ACTION printouts may differ from these figures because of 
modifications to the program. The version used to generate these figures 
operated under the CDC 6000 computing system, used the BCD character set 
and printed from single-precision floating point data. 

Problem Definition Printout 

Figures 13 through 21 present illustrations of typical printout for 
the problem description. This description is the program's interpretation 
of the input data. Figure 13 illustrates how the program title and 
control parameters are displayed. Figures 14 through 20 present printouts 
of input entered via data blocks. Included in the problem definition 
printout is a summary of descriptive parameters and storage locations. 

Figure 21 illustrates this printout which concludes the problem definition 
output. This summary terminates ACTION printed output whenever errors 
are detected in input data. 

Printout of Response Data (Equilibrium Positions) 

ACTION provides a time history of the response of a structure subjected 
to various loads and initial conditions. This history is in the form of 
printed records of equilibrium configurations. It includes a description 
of solution progress and accuracy, joint motions, element and system 
energies, and, as an option, stresses. Printouts of configurations are 
provided at each load step or each time step needed for numerical 
integration of the equations of motion, provided IPRT equals zero or one. 

A record of progress made in minimizing the energy function for each 
time increment precedes any printout of joint and element information. 
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Figure 13. Printout of Problem Title and Control Parameters 
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Figure 14. Printout of Joint Coordinate Data 
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Figure 15. Printout of Material Data 
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Figure 16. Printout of Line Element Data 
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Mass Data 
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Figure 18. Printout of Initial Conditions 
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Figure 19. Printout of Function Tables 
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Figure 20. Printout of Forcing Funtions 
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Figure 21. Printout of Summary of Parameters and Storage Location 
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This record provides data indicating the path taken in arriving at the 
solution of the equations of motion of the system. The IINK parameter 
of the first Control Parameter Card controls the amount of minimization 
data provided. Figures 22 through 24 illustrate printout for IINK values 
of 0, 1, 2, and 3, respectively. 

Information regarding solution accuracy, the state of joints, and 
stored and dissipated energies is always printed at the end of each time 
step. Figure 25 and 26 exhibit examples of this printout. 

A measure of solution accuracy is provided by the zero, half, and 
full time error data of Figure 25. This data refers to a normalized 
equilibrium imbalance obtained for each degree of freedom at the 
beginning, middle and end of the time step, respectively. The equilibrium 
imbalance is taken to be the variation of the energy function with 
respect to a small change in the displacement component, multiplied by 
the displacement component and divided by the current value of the 
energy function. The error measures displayed at the top of Figure 25 
are the maximum absolute values of the individual degree of freedom 
errors. From these a single measure of error is computed as the difference 
between the half time error and the average of the zero and full time 
errors. (If this value is larger than an upperbound specified on the 
second Control Parameter Card the time step is reduced by 0.5 and the 
integration for the time step is repeated. If the error is smaller than 
a specified lower-bound the time step is increased by 1.50 for the next 
integration cycle. In addition the time step is not permitted to be less 
than a specified minimum and the time step will also be reduced if the 
full time error is greater than the upper-bound requirement.) 
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Figure 22. Printout of Minimization Progress (IINK = 0) 
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Figure 23. Printout of Minimization Progress (I INK = 1) 
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Figure 25. Printout of State of Joints 


Figure 25 also displays printout describing the state of joints for 
the time indicated at the top of the page. The current displacement, 
velocity, acceleration, jerk and applied forcing function is defined 
for each joint degree of freedom. DX, DY, and DZ refer to displacement 
degrees of freedom along the global, X, Y, and Z axes, respectively; 

RX, RY, and RZ refer to rotational degrees of freedom about the global 
X, Y, and Z axes, respectively. 

The total response of the structural system is summarized by the 
system energy evaluations printed at the top of Figure 25, The total 
input energy of the system is the total amount of work done on the system 
since the integration began. This total energy is then broken down 
first into the current system energy which is the sum of the current 
elastic strain energy and kinetic energy, and then into simply the current 
elastic strain energy. Finally, the current dissipation rate is printed. 
This is the rate at which strain energy was lost due to yielding during 
the current time step. 

Figure 26 shows the output defining how the strain energy stored 
as elastic energy and lost due to yielding during the current time 
step is attributed to individual elements. By examining the percentage 
of the total system strain energy contributed by each element, the 
contribution of each element to the total structural stiffness becomes 
readily apparent. Plus and minus errors on this printout provide 
measures of the accuracy of numerical integration over the element 
cross section. Errors are reported in percent. 

Stress and strain-rate printout is provided when the ISINK parameter 
of the first Control Parameter Card is set equal to 1. When ISINK = 0, 
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Figure 26. Printout of Element Energy Status 


this printout is suppressed. When ISINK = 2, membrane element generalized 
force printout is produced as well as all printout available at the ISINK = 1 
level. Figure 27 presents a typical printout. Stress and strain-rate 
values are provided for up to 16 locations throughout a cross-section. 

These locations are referred to as stress points and are illustrated in 
Figures 7 through 11. 

Pages providing stress output also include failure checks for 
fracture, buckling, deformation, dissipated energy and shear yield. 

ACTION simulates fracture at individual quadrature points if the axial stress 
at a point would, exceed the ultimate stress of the material. Fracture 
is indicated by "1.0" in the "FAILURE" columns of the printout. The 
buckling check flags compression forces exceeding eight-tenths of the 
axial buckling force as calculated from element input data. The buckling 
force is included in the printout of member section properties. The 
deformation check flags relative rotations within elements exceeding 0.3 
radians. The dissipated energy check warns of the existence of an energy 
distribution not properly accounted for in the formulation of the 
energy function. The deformation model neglects shear deformations in 
the plastic region. The check on shear yield warns the user when this 
approximation may lead to significant errors. These latter four checks 
are printed following the printout for individual stress points.. If 
ISINK > 0, the actual values used in the checks are printed. Failure 
will then be indicated by M 1 M In the appropriate failure column. If 
ISINK = 0, the checks are printed only when the analysis determines 
failure.. These checks are not fatal but merely serve as a warning to 
the user. 
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Figure 27. Printout of Stress Psesults 


Error Messages 


Diagnostic error messages describing input errors and solution 
errors are printed as errors are encountered during execution of the 
program. These error messages are self-explanatory. 
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SAMPLE PROBLEM 


This section, provides additional information on the use of scratch 
files, restart tapes, etc. in attempting to run a sample problem using 
a CDC procedure file and/or an IBM JCL file. The objective in providing 
such sample input and selected output, as in Figures 25 through 27, is 
to provide the users with some means to validate their version. The 
selection of the sample problem was motivated by the availability of an 
independent solution for correlation purposes. Figure 28 provides the 
input files for the sample problem. 

A discussion of this and a few other problems that demonstrate the 
overall capability of the ACTION program may be found in refs. [5] and 
[ 6 ]. 
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1. . Instructions for Input: 


1.1 CDC Machine: 


The following is a procedure file that may be used to excute a 


data deck read from cards on a CDC machine. 


JOB CARD 
USER CARD 
CHARGE CARD 
COPYCR(DATA) 

REWIND (DATA) 

GET ( BACT /NA) 

GET (TAPE 10) 

REWIND (BACT , TAPE10) 

GET (NMACFTN/UN=LIBRARY ) 

LDSET (PRESETA=NGINF , LIB=NMACFTN) 
BACT (DATA) 

REPLACE (TAPE10) 

GOTO , 1 . 

1, RETURN (BACT) 

RETURN(TAPEIO) 

REWIND (DATA) 

COPYSBFDATA, OUTPUT) 

/EOR 


DATA DECK (see next page) 


/EOF 


Note that BACT is a binary file of the ACTION source which is presumably 
residing on user's permanent disk space. The cards involving TAPE10 are 
required only for modifying a restart tape. 
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Figure 28. Input Data for an Impulsively Loaded Clamped Beam 
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1.2 IBM Machine: 


The following is a Job Control Language file that may be used to 
execute a data deck read from cards on an IBM machine. 


JOB CARD 
PRIORITY CARD 
/*JOBPARM LINES=XX 

/ /STEP1 EXEC PGM=LMACT , TIME=XX , REGION= 1000K 
//STEPLIB DD DSN=AXXXXX. ACTION, DISP=SHR 
//GO.FT12F001 DD DSN=&TEMP,UNIT=SYSDA, 

// SPACE=(TRK, (30,10)) ,DCB=(RECFM=VSB,LRECL=8000,BLKSIZE=8004) 
//GO.FT10F001 DD DSN=AXXXXX . ACTIO , UNIT=SY SDA , 

// SPACE= (TRK , (60,60)) ,VOL=SER=XXXXX,DISP= (OLD, KEEP) , 

// DCB=(RECFM=VSB ,LRECL=8000 ,BLKSIZE=8004) 

//FT06F001 DD SYSOUT=A 
//FT05F001 DD * 


DATA DECK (see previous page) 


/* 

// 


Note that LMACT is the load module member name in the data set AXXXXX. ACTION 
previously created. The GO.FT10F001 card is required only in the event that 
a restart tape, is being written. 


2. Sample Problem Details : 


Problem Title: Impulsively Loaded Clamped Beam 

Comments: This problem uses the LINE Element IE to model a 

beam with a solid rectangular cross-section. Since the beam under- 
goes bending about only one of its principal axis this element is 
preferred over the other LINE element SORE which should be used 
for the more general case of bending about both prinicpal axes. 

Figure 29 compares the predicted response using two versions 
of the Newmark-Beta method with that known for experiments [7]. 
Results indicate sensitivity of the response to assumed values 
of the 6 and y parameters of the Newmark-Beta Method. 
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Figure 29. Displacement Response of the Impulsively Loaded Clamped Beam 
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